Vamos a ver si el que el inmueble esté cerca del metro tiene algún efecto en el precio. Para ello nos bajamos la geolocalización de las paradas de metro de Nomecalles
bocas <- readOGR(dsn="./data/Transportes y comunicaciones_ Metro (bocas)/",layer="bocas")
## OGR data source with driver: ESRI Shapefile
## Source: "./data/Transportes y comunicaciones_ Metro (bocas)/", layer: "bocas"
## with 740 features
## It has 7 fields
## Integer64 fields read as strings: UTM_Y
Y las transformamos a lat/long
WGS84 <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
bocas2 <- spTransform(bocas,CRSobj=WGS84)
Ponemos la latitud y la longitud en los datos
cc <- coordinates(bocas2)
bocas2@data$lon <- cc[,1]
bocas2@data$lat <- cc[,2]